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Abstract. Verification of theoretical predictions of an oscillating behavior of the fine-structure constant a (= /he) 
with cosmic time requires high precision Aa/a measurements at individual redshifts, while in earlier studies the 
mean Aa/a values averaged over wide redshift intervals were usually reported. This requirement can be met 
via the single ion differential a measurement (SID AM) procedure proposed in Levshakov et al. (2005). We apply 
the SIDAM to the Fell lines associated with the damped Lya system observed at Za_hs~ 1.15 in the spectrum 
of HE 0515-4414. The weighted mean {Aa/a) calculated on base of carefully selected 34 Fell pairs {A1608, X} 
{X = 2344,2374, and 2586 A) is (Aa/a) = (-0.07 ± 0.84) x lO"*^ (la C.L.). The precision of this estimate 
improves by a factor 2 the previous one reported for the same system by Quast et al. (2004). The obtained result 
represents an absolute improvement with respect to what has been done in the measurements of Aa/a . 
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1. Introduction 

One of the most important physical constants, the 
Sommerfeld fine-structure constant a, is subjected to more 
and more precise measurements in both modern labora- 
tory experiments and astronomical observations due to 
fundamental role which a plays in quantum electrody- 
namics effects and in electromagnetic and optical proper- 
ties of atoms. An unprecedented accuracy of ~ 0.1 x 10~^ 
(the relative error) has been recently achieved in the lab- 
oratory measurements of the atomic transitions in alkali 
atoms (Banerjee et al. 2004). This result, combined with 
sub-ppb (parts per billion) measurements of the atomic 
masses (Bradley et al. 1999), enables a high-precision de- 
termination of the fine-structure constant at the present 
time, i.e. at the cosmological epoch z = 0. 

The question whether or not the value of a var- 
ied at different cosmological epochs can be answered 
only through spectral studies of extragalactic objects. 
Theoretically the effects of inhomogeneous space and time 
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* Based on observations performed at the VLT Kueyen tele- 
scope (ESO, Paranal, Chile). The data are retrieved from the 
ESO/ST-ECF Science Archive Facility 



evolution of a were considered by Marciano (1984) and 
Mota & Barrow (2004). Most recently Fujii & Mizuno 
(2005) and Fujii (2005) suggested a damped-oscillation- 
like behavior of a as a function of cosmic time t. It is ap- 
parent that to study such irregular changes in a, we need 
to achieve high precision in the measurements of Aa/a at 
individual redshifts, contrary to the averaging procedure 
over many redshifts which is usually used to decrease un- 
certainties of the mean values (Aa/a) (see, e.g.. Murphy 
et al. 2004, and references therein). 

Contemporary astronomical observations still cannot 
provide accuracy comparable to the laboratory results. Up 
to now, the errors of the fractional deviation of the fine- 
structure constant, Aa/a = (a^— ao)/ao, recovered from 
astronomical spectra are at the level of 10~^ (or parts per 
million, ppm)^. For instance, individual Aa/a values ob- 
tained by Chand et al. (2004) are known with the accuracy 
of a few ppm. The uncertainty of the same order of magni- 
tude, cr^Aa/a) — 1.7-1.9 ppm, was achieved in our previ- 
ous study of the Zahs = 1-15 system towards HE 0515-4414 
(Quast et al. 2004, hereafter QRL). In both cases the stan- 
dard many-multiplet (MM) method (Webb et al. 1999; 



Here ao and denote the values of the fine-structure con- 
stant at epoch 2 = and in a given absorption (or emission) 
line system at redshift z, respectively. 
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Table 1. ESO UVES archive data on the quasar HE 0515-4414 



Exp. 


Set- 


Slit 




Time, 




*~v 001 T1 rr 




P, mb 




T-'ttmtt'ci m m o 


No. 


tine 


arcsec 


v-m-d 


li:ni 


h:m 


arcsec 






K 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


lb 


346 


0.8 


1999-12-14 


05:27 


1:15 


0.46-0.56 


11.8-12.8 


c\r\r\ A r\c\-\ r\ 

990.4-991.0 


0.2 


dO.A-9022 


1 r 


OoU 


n 7 
u. * 


















2b 


346 


0.8 


1999-12-14 


06:43 


1:15 


0.56-1.04 


11.5-11.9 


c\c\r\ -1 r\c\r\ o 

990.1-990.3 


0.3 


dO.A-9022 




Pi Sin 

OoU 


n 7 
u. * 


















3b 


346 


0.8 


1 c\c\c\ 1 r\ 1 r" 

1999-12-15 


04:47 


1:30 


0.51-1.84 


10.1-11.9 


c\c\(\ r\ nnr» o 

990.0-990.8 


0.006 


dO.A-9022 


Ol 


OOU 


n 7 


















4b 


346 


0.8 


C\ l\ f\ l\ -1-1 -1 ^ 

2000-11-17 


06:42 


1:15 


0.34-0.53 


13.5-13.8 


990.1-990.3 


1.0 


DD.A-0212 


Arc 


OOU 


u.o 


















5b 


346 


0.8 


2000-11-18 


07:25 


1:15 


0.55-0.85 


13.2-13.7 


989.1-989.2 


0.8 


66.A-0212 


5r 


580 


0.8 




07:25 














6b 


346 


0.8 


2000-11-19 


07:27 


1:15 


0.58-0.80 


13.2-14.0 


989.2-989.5 


0.2 


66.A-0212 


6r 


580 


0.8 




07:27 














7b 


346 


0.8 


2000-12-25 


06:13 


1:15 


0.55-0.88 


15.8-16.5 


989.9-990.0 


0.5 


66.A-0212 


7r 


580 


0.8 




06:13 















Seeing, S, temperature, T, and air pressure, P, limiting values listed, respectively, in Cols. 7, 8 and 9 are taken from 
the ESO/ST-ECF Science Archive Ambient Conditions Database at ,http:/ /archive.eso.org/| 



Dzuba et al. 1999, 2002) has been used. Further modifi- 
cation of the MM method (Levshakov 2004; Levshakov et 
al. 2005, hereafter LCMD) resulted in a new methodology 
for probing the cosmological variability of a on base of 
pairs of Fell lines observed in individual exposures from 
a high resolution spectrograph (henceforth referred to as 
SIDAM - single ion differential a measurement). The ba- 
sic idea behind SIDAM was to avoid the influence of small 
spectral shifts due to ionization inhomogeneities within 
the absorbers and due to non-zero offsets between differ- 
ent exposures. The individual offsets can affect the shape 
of the line profiles during rebinning and coadding pro- 
cedures which are usually applied to combine exposures 
together to increase signal-to- noise, S/N, ratio (examples 
are given in LCMD). 

In the present paper we show that the SIDAM can 
provide a sub-ppm precision in a single redshift Aa/a 
measurement and that this level of accuracy is caused by 
intrinsic instrumental imperfections and systematic errors 
inherited from the uncertainties of the wavelength scale 
calibration. 

While working at sub-ppm level, we confront the 
problem of sub-pixel centering. As shown by David & 
Verschueren (1995), this problem can be properly treated 
only in two cases: either (i) the exact analytical func- 
tion describing the observed line profile is known a priori, 
or (m) the observed line profile is intrinsically symmetric. 
Neither of these conditions is fulfilled for metal absorption 
lines observed in QSO spectra. This means that in general 
the line centering must be handled with special care. We 
discuss this problem in detail in Sects. 2 and 3. The results 
of the SIDAM are presented in Sect. 4. Our conclusions 
are given in Sect. 5. 



Table 2. ESO UVES archive data on ThAr lamp used to 
calibrate spectra of HE 0515-4414 



ThAr 


Date, 


Time, 


T,°C 


P, mb 


No. 


y-m-d 


h:m 






1 


1999-12-14 


08:17 


11.8 


990.0 


4 


2000-11-17 


07:58 


13.7 


990.1 


5 


2000-11-18 


08:41 


13.3 


989.1 


6 


2000-11-19 


08:44 


14.1 


989.3 


7 


2000-12-25 


07:29 


16.5 


990.0 



Temperature and air pressure are taken from the 
ESO/ST-ECF Science Archive Ambient Conditions 
Database at .http: / /archive. eso.org/. 



2. Observations and Data Reduction 

We analyze high quality spectra of the bright intermedi- 
ate redshift quasar HE 0515-4414 (zom = 1-73, B = 15.0; 
Reimers et al. 1998). The observations were acquired with 
the UV- Visual Echelle Spectrograph (UVES) at the VLT 
8.2 m telescope at Paranal, Chile, and the spectral data 
were retrieved from the ESO archive. The selected expo- 
sures are listed in Table 1 . The spectra were recorded with 
a dichroic filter which allows to work with the blue and red 
UVES arms simultaneously as with two independent spec- 
trographs, and the CCDs were read-out in 1 x 1 binned 
pixels (spatial X dispersion direction). The standard set- 
tings at central wavelengths A346 nm and A580 nm were 
used for the blue and red arms respectively (marked by 
symbols 'b' and 'r' in Col.l of Table 1 (along with the ex- 
posure number). From the blue spectra we used only order 
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Table 3. Original pixel sizes along the selected echelle 
orders 



-100 



-150 



19b 8r 9r lOr 

echelle order 



19r 



20r 



Fig. 1. Doppler shifts of two calibrated ThAr lines brack- 
eting the positions of the Fe ii absorption complex (Fig. 2) 
for exposure No.l (Table 2). Echelle orders (Table 3) are 
indicated on the bottom axis ('b' and 'r' denote the blue 
and red arms, respectively). The ThAr lines used are 
the following (wavelengths in A): 3457.070 & 3462.851 
(19b), 5039.230 & 5041.122 (8r), 5101.129 & 5111.278 (9r), 
5122.499 & 5125.489 (lOr), 5559.891 & 5564.201 (19r), and 
5595.063 & 5599.654 (20r). The peak-to-peak variations of 
A^^noise ~ 100 m s~^ Set a lower limit to the error of the 
Fell relative position measurements 



19, and from the red spectra orders 8, 9, 10, 19 and 20^, 
where Feii lines suitable for the Aa/a measurement are 
observed. All selected Fell lines are located close to the 
central regions of the corresponding echelle orders. This 
minimizes possible distortions of the line profiles caused by 
the decreasing spectral sensitivity at the edges of echelle 
orders (LCMD). 

Other details of the observations are also presented in 
Table 1. Cols. 2 and 3 give setting and slit width used. 
Cols. 4 and 5 list the date and the starting time of the 
exposure whose duration is indicated in Col. 6. The am- 
bient conditions are characterized by minimum and max- 
imum values of the seeing, 5, temperature, T, and air 
pressure, P (Cols. 7-9). During an integration these quan- 
tities were varying between the indicated limits. Col. 11 
shows that the data were obtained for two observational 
programs, 60.A-9022 and 66.A-0212. To characterize the 
quality of the obtained spectra, we introduce a quality fac- 
tor, K (Col. 10), defined as k = 1/(A>S • AT ■ Apj and 



^ These are the sequential numbers of the echelle orders 
which are used throughout the paper. The correspondence be- 
tween them and echelle orders is as follows: 19b=135, 8r=121, 
9r=120, 10r=119, 19r=110, and 20r=109. 



Echelle 


AA (A) 


A(Fc 


II) complex 


(A) 




l-'ixd .siz(> (luA) 


I'ixo 


. size (km s 




HliiP firm 


'•!44q 'Mfiq 




3459 




#19 


19.5-16.0 




1.55 




vpH arm 


5030-5050 




5040 




#8 


22.7-20.5 




1.29 




■ppfl arm 


5095-5115 




5104 




#9 


20.4-17.9 




1.13 




red arm 


5113-5133 




5123 




#10 


23.2-21.0 




1.30 




red arm 


5552-5572 




5562 




#19 


22.9-20.5 




1.17 




red arm 


5581 5601 




5591 




#20 


25.3-23.2 




1.30 





expressed in fractions of the 4th exposure which has the 
highest quality. The selection of variables A5, AT, and 
AT is heuristic in some sense but it is relevant to the 
problem of the line centering. For instance, a change of 
AP = 1 millibar (or a change of AT = 0.3°C) induces an 
error in radial velocities of ~ 50 m s~^ (Kaufer, D'Odorico 
& Kaper 2004). Problems may also occur because of seeing 
variations which change the intensity of the QSO signal 
during an integration, and, hence, the instabilities of the 
spectrograph may be sampled in different ways. Other ef- 
fects like vibration, thermal drift, changes of the grating 
spacing due to temperature fluctuations, unequal illumi- 
nation of the spectrograph grating and collimator by the 
laboratory reference source and starlight, varying width 
of the instrumental profile, and scattered light may pro- 
duce additional Doppler noise at the level of a few m s~^ 
(Griffen & Griffen 1973; Brown 1990; Meyer 1990; Gulliver 
et al. 1996). These deleterious effects do not permit very 
precise measurements of the line positions even from ex- 
tremely high S/N spectra. However, these problems can 
be considerably reduced using fiber fed spectrographs like, 
e.g., a configuration of UVES-I-FLAMES which provides 
the radial velocity precision better than 100 m s^^ for 
the stars with V magnitudes in the range 14-18 (Melo et 
al. 2004; Bouchy et al. 2005). 

We used the UVES pipeline (the routines implemented 
in MIDAS-ESO data reduction package for UVES data) to 
perform the bias correction, inter-order background sub- 
traction, flat-fielding, correction of cosmic rays impacts, 
sky subtraction, extraction of the orders in the pixel space 
and wavelength calibration. A modified version of the rou- 
tine 'echred' of the context ECHELLE inside MIDAS was 
used to calibrate in wavelength the echelle spectra without 
rebinning. In this way we have the reduced spectra with 
their original pixel size in wavelength. 

The residuals of the calibrations are rather small 



amounting to CTr 



< 



1 mA. The observed wavelength 



4 
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Table 4. The instrumental profile widths (FWHM) in km s ^ and the mean signal-to-noise ratios per pixel (given in 
parentheses) at the continuum level at the positions of Fe ii lines in scientific exposures 



Echelle 
order 



#1 



#2 



Scientific exposures 
#3 #4 #5 



#6 



#7 



19b 
8r 
9r 
lOr 
19r 
2l)r 



5.60 



+°o:iS(3i) 



5.60l°:l°(31) 



5+0.07 
'-0.07 

fO.05 

-0 



(46) 



S+0.07 



(50) 



4.77lHi:(51) 4.77«:;»(57) 



4.841^^9(47) 
4.78l°;°3(47) 



".•■!) 



4.84«;«^(43) 
4.78l°:°l(54) 

i.8:r"-;:^(i!)) 



5.6oi°:iS(27) 

4.88lH^(43) 

4.77i;::;«(49) 

4.84: 



S:iJ^9(42) 



4.781°;°! (49) 



i.,s:i-t 



15) 



5.60l°:l°(31) 

5.50l°:l°(56) 



+0.09 
0.09 



5.36 



5.46«;l2(49) 
5.391° ^^(53) 



-0.11 

■■n.o: 



(ii) 



5.60l°:l°(31) 



5.40 



+0.10 



(54) 



5.60l°:l°(28) 

5.50l°:^S(51) 



(44) 5.43i;;:»;:(47) 5.47l°:««(52) 
5.46l«:f2(52) 5.54l°;}l(50) 



5.45l°;i(54) 

"..17-' " I! (GO) 



5 fil+O-27 
^•0^-0.27 



(59) 

--0.11 



.11 ^ 



5.601°:!° (35) 



5.50 



+0.06 



(51) 
(45) 



5.47l,,oB 
5.461°:°^ (53) 
5.45l°;ll(54) 

5. 12-'-;;;- (18) 



scale of each spectrum was transformed into vacuum, he- 
liocentric wavelength scale (Edlen 1966). 

Following the same procedure as for the QSO expo- 
sures, we calibrated in wavelength the spectra of the ThAr 
arcs. This allowed us to estimate the 'Dopplcr noise' by 
measuring the random velocity shifts of the ThAr emis- 
sion lines. For each order we selected two well-exposed 
ThAr emissions which bracket the positions of the Fell 
absorption complex. In Fig. 1 we show the Dopplcr shifts, 
Atinoise, mcasurcd with the calibrated ThAr lamp No.l 
(Table 2). For both lines of the arc the peak-to-peak vari- 
ations of Az^noise are about 100 m s~^ . This gives us a 
lower limit to the error of Fell positional estimations. 

Usually when a precise absolute wavelength scale is 
required, calibration exposures are taken before and af- 
ter the scientific exposures. As discussed above, varying 
ambient weather conditions may introduce different ve- 
locity offsets in the lamp and QSO spectra if they were 
not obtained closely in time. However, this uncertainty 
affects only an absolute calibration whereas in Aa/a es- 
timations we are dealing with differential measurements. 
In fact, we measure the relative Fell positions with re- 
spect to the Fell A1608 hnc (LCMD). Since the SID AM 
uses iron lines obtained simultaneously in one exposure, 
a systematic offset caused by shortcomings of the calibra- 
tion procedure should be canceled out for the blue and red 
arms spectra. This is not the case if unexpected mechani- 
cal instabilities occur during integration (see Sect. 3 below 
and an example in LCMD where effects of the mechanical 
instabilities are discussed). For our data, the first three 
QSO exposures were calibrated with the same ThAr lamp 
No.l (Table 2). The 4th to 7th QSO exposures were cali- 
brated with the corresponding ThAr lamps No. 4-7 taken 
immediately after the scientific exposures (cf. Tables 1 and 
2). 

Our further concern was nonlinearity of the wavelength 
scale. Off-plane design in echelle spectrographs like UVES 
is introducing a pronounced curvature of the orders (e.g. 
Ballester & Rosa 1997). As a result, the wavelength scale 
is not linear in the sense that the original pixel width in 
wavelength decreases with increasing pixel number along 
the frame. Table 3 illustrates this behavior for the selected 
echelle orders. In Col. 2, upper row for each echelle order 



indicates the wavelength range which includes the Fe ii 
absorptions, whereas the corresponding pixel sizes at the 
starting and end points of this range are given below. Col. 3 
shows the pixel size at the center position of the Fe ii com- 
plex (lower and upper rows, respectively). The Fell com- 
plex is depicted in Fig. 2. 

The difference between pixel sizes at the edges of the 
wavelength ranges shown in Fig. 2 and Table 3 is about 
18% for order 19 (blue arm) and 8-12% for the red arm 
orders. This introduces an artificial inclination of the local 
continuum level which may affect the relative positions of 
the Fe ii lines at the sub-pixel level. Therefore we normal- 
ized the registered photocounts per pixel to the original 
pixel size in wavelength and determined the local continua 
for each Fe ii complex following linear regression procedure 
from LCMD. The uncertainties of the calculated continua 
are less than 1% in all Fell regions. An example of the 
continuum fitting is shown in Fig. 2 for the 4th QSO ex- 
posure. 

Finally, we measured the FWHM of the instrumental 
profile for each individual echelle order. The instrumental 
profile is dominated in our case by the slit width which 
was different for the blue and red arms in the first three 
exposures (Table 1). The FWHM values calculated from 
the narrow lines of the arc spectra are given in Table 4. 
Our further results were obtained with the mean FWHM 
and signal-to-noise ratio values listed in Table 4. 

3. Concordance of the Feii profiles 

The spectra of HE 0515-4414 reveal a multi-component 
complex of metal absorption lines associated with a sub- 
damped Lya (sub-DLA) system at Zabs =1.15 (dc la Varga 
et al. 2000; Quast et al. 2002; Reimers et al. 2003). The 
radial velocities of the absorption components of the Fe ii 
complex span 660 km s~^ (QRL, Table 2). Not all of 
them are seen, however, in the individual exposures, the 
weakest ones were detected only in the co-added spectra 
(QRL, Fig. 1). Therefore, for our purpose we chose two 
sub-systems at -Zabs= 1.150965 and 1.149092 which ex- 
hibit the most pronounced absorption features and, thus, 
provide the most accurate determination of their relative 
positions from individual exposures. The atomic data used 
for the analysis are given in Table 5. 
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5 



blue arm #19 



Fell 1608 



in 



X 



3449 3450 3451 3452 3453 3454 3455 3456 3457 345S 3459 3460 3461 3462 3463 3464 3465 3466 3467 3468 3469 



— . .... ....... ,,. .1.' -11 Ml -I'^.i ,M n an. 1,.- r .r.'-i'' -'-,11 












red arm #8 






Fell 2344 











> 1 < 


1 . . . . 



5030 5031 5032 5033 5034 5035 5036 5037 5038 5039 5040 5041 5042 5043 5044 5045 5046 5047 5048 5049 5050 



red arm #9 



Fell 2374 





5095 5096 5097 5098 5099 5100 5101 5102 5103 5104 5105 5106 5107 5108 5109 5110 5111 5112 5113 5114 5115 



' ' I ' ' ' ' I ' ' ' ' I 



red arm #10 



Fell 2382 



5113 5114 5115 5116 5117 5118 5119 5120 5121 5122 5123 5124 5125 5126 5127 5128 5129 5130 5131 5132 5133 
' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 



red arm #19 
"i^^i 



Fell 25B6 




5552 5553 5554 5555 5556 5557 5558 5559 5560 5561 5562 5563 5564 5565 5566 5567 5568 5569 5570 5571 5572 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 










red arm #20 






Fell 2600 








1 i' '..'1 i . "i'-j- 






5581 5582 5583 5584 5585 5586 5587 5588 5589 5590 5591 5592 5593 5594 5595 5596 5597 5598 5599 5600 5601 

Wavelength (A) 

Fig. 2. Unnormalized portions of tlie HE 0515-4414 spectra (intensities are in arbitrary units) with the Fell lines 
(2abs= 1.15) obtained from the 4th exposure (Table 1). All Fell lines are located near the center of the orders, 
(indicated at the left side of each panel). Dots with error bars arc the mean intensities and their la cirrors (calculated 
in the ranges marked with horizontal lines) used to estimate the local continua by means of the linear regression 
analysis (LCMD). The uncertainty of each local continuum level is less than 1% 



In both sub-systems there arc no isolated and symmet- 
ric absorption components which could be used for the 
line centering. But, as explained in LCMD, we are mainly 
interested in the position of an absorption complex as a 
whole to calculate its radial velocity shift with respect to 
the reference line Fell A1608 — the only line among Fell 
transitions which has a negative sensitivity coefficient Q 
(Table 5). 

In general, the position of the line blend is sensitive 
to the relative strengths of the partly resolved and un- 
resolved (hidden) components as well as to the shape of 
the instrumental profile. The effect is most noticeable in 



the ease of optically thick and narrow hidden components 
due to their different saturation depending on the line 
strength (Levshakov 1994; Levshakov & D'Odorico 1995; 
QRL). As seen from Table 5, the strengths of the iron 
lines, log/Avacj most frequently observed in QSO spectra 
(Fell AA2374,2382) differ by an order of magnitude. 

To handle the blend we have to adopt a specific model. 
Since we are dealing with only one ion, all Fell lines 
from an intervening absorber must have similar profiles. 
A model not matching all observational profiles simulta- 
neously either indicates blended lines from different inter- 
vening systems or that some data points are outliers. This 
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Table 5. Atomic data of the Feii transitions", and the 
sensitivity coefBcients Q*". Estimated errors are given in 
parentheses 



500 



iVli L . 


\ 

/*vac7 


J 




iOg / '^vac 


No.= 


A 








lu 


2600.1725(1) 


0.23878 


0.035(4) 


2.79 


lu 


2586.6496(1) 


0.06918 


0.039(4) 


2.25 


2u 


2382.7642(1) 


0.320 


0.035(4) 


2.88 


2u 


2374.4603(1) 


0.0313 


0.038(4) 


1.87 


3u 


2344.2130(1) 


0.114 


0.028(4) 


2.43 


8u 


1608.45080(8) 


0.0580 


-0.021(5) 


1.97 



"based on the compilation of Murphy et al. (2003); 
''defined in Sect. 4; '^multiplet numbers from Morton (2003) 



allows us to check concordance of the Fell data before 
evaluating Aa/a . 

To define the model wc will assume, following LCMD, 
that: (1) the number of subcomponents ng is fixed for all 
Fe II lines under study; (2) the Doppler bi parameters are 
identical for the same ith subcomponents; (3) the rela- 
tive intensities of the subcomponents rij and (4) the rela- 
tive radial velocities 5vij between the subcomponents are 
fixed. 

Then, an Fe ii blend can be described by the sum of rig 
Voigt functions: 



(1) 



i=l 



where ri is the optical depth at radial velocity v within 
the £th Fell line (£ = 1,2,..., L), Ni is the column density 
of the main component, r^.i — Ni/Ni, vg is the position of 
the main component in the line £, and dvi^i = 0. Here L 
is the total number of Fell lines involved in the analysis. 

The model is fully defined by specifying p = 3ns + L — 
1 parameters: A^i, {bijZi, {^Vi,i}Zi\ {ri,i}Zi\ and 
{ve}f^i. All these parameters are components of the pa- 
rameter vector 9 = {6i, 62, ■ ■ ■ , Op}. To estimate 6 from the 
combined Fe 11 profiles, we minimize the objective function 



-J L mi 



(2) 



=1 j=i 



where J^f^^ is the observed normalized intensity of the 
spectral line £, aij is the experimental error within the jth 
pixel of the line profile, and p = J2e ~ P = M — p is 
the number of degrees of freedom. Tlj{0) is the calculated 
intensity convolved with the corresponding spectrograph 
point-spread function (Table 4). 

The total number of data points involved in the anal- 
ysis is M = 2704 (zabs = 1.150965) and M = 428 ( Zabs = 
1.149092). All selected Fell profiles are free from cosmic 
rays and telluric absorptions. In accord with QRL, we took 
TT-s = 8 (p = 65) and = 2 [p = 47) to define our models 
for the former and the latter sub-systems. We found that 
the models chosen described both sets of data adequately: 
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Fig. 3. Same as Fig. 1 but for the Fell lines from the 
Zahs = 1.150965 system. The exposure numbering is given 
in accord with Table 1, whereas the line strengths and 
wavelengths - according to Table 5. The relative velocity 
offsets are based on the 8-component model estimations 
(see text for details) 



the xt per degree of freedom was about 1 in both cases. 
Thus, each ensemble of Fell lines is self-consistent, free 
from outliers and blending with other lines. 

This preliminary analysis revealed, however, an unex- 
pectedly large shift between the blue and red arms for the 
3rd exposure (short-dashed line in Fig. 3). It might be a 
chance coincidence, but the spectral data obtained from 
the 3rd exposure have the lowest quality factor k = 0.006 
as compared with unity (see Table 1). The quality factors 
for other exposures are much larger and do not exhibit 
significant variations in their values. Similar offsets be- 
tween the blue and red arms were noted for the UVES 
spectra of Q 1101-264 in LCMD where we ascribed them 
to mechanical instabilities. Such exposures were excluded 
from further Aa/a measurements since they may induce 
a mock signal. 

Fig. 3 shows one more problem: there are too large 
peak-to-peak variations of the relative velocity offsets, Au, 
and we note a clear dependence of At; on the line strength. 
This behavior can be corrected by excluding saturated 
lines (like Fe 11 A2382 and A2600) from the analysis, and by 
increasing the number of sub-components to account for 
the influence of the unresolved (hidden) blends on the line 
centering. The number of sub-components is controlled by 
two factors: (i) the value of Xmin should be close to unity, 
while (ii) the dispersion of Au should decrease but not 
differ much from that deduced from ThAr lines shown in 
Fig. 1. Fig. 4 (upper panel) shows the results obtained for 
a model with rig = 13 sub-components [p = 62). 
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Fig. 5. Individual exposures (labeled as 'e#') of the Fe ii lines selected from the spectrum of HE 0515-4414 to estimate 

Aa/a. Normalized intensities arc indicated by histograms. The over-potted synthetic profiles (smooth curves) arc 
calculated from the joint analysis of all iron lines. The zero radial velocity is fixed at z = 1.150965. The minimization 
procedure gives Xmin = 0-9 P^r degree of freedom {u = 2642). The dashed vertical lines mark positions of the sub- 
components from Table 6 



The estimation of the uncertainties of the best-fitting 
model parameters in case of many-dimensional parame- 
ter space requires a special approach. When p ^ 1, the 
parameters are as a rule correlated, the parameter space 
near the minimum of the objective function has a com- 
plicated topology and the Ax^ levels for the confidence 
regions cannot be accurately defined. In such cases Monte 
Carlo simulations are recommended instead of the for- 
mal inversion of the covariancc matrix (sec, e.g.. Press et 
al. 1989, Chapter 15). We used the bootstrapping resid- 
uals method to estimate the standard error (arms) of the 



best-fitting parameters. The trial datasets are created in 
the following way. Suppose 6 is the estimated parameter 

vector and rej = — are the residuals. Then a 

new dataset J-^j^' is made by sampling independently 
from re J for each line £, yielding r|^^-, and by setting 
jr°bs* ^ J^^^l{0) + r} .. New values 6* of the parameter 
vector are then computed from the bootstrap data in the 
same way that 9 was computed from the original data, 
i.e., by least squares. 
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Fig. 6. Same as Fig. 5 but for the sub-system at -2abs= 1.149092. The minimization procedure gives Xmin = 0-^ P^^ 
degree of freedom {v = 399). The dashed vertical lines mark positions of the sub-components from Table 7 



The model parameters and their uncertainties com- 
puted on the base of 30 bootstrap samples are given in 
Tabic 6. The pcak-to-pcak variations arc now diminished 
by a factor of 2-3 and are equal to 200 m s^^ in average. 

The lower panel in Fig. 4 shows the peak-to-peak vari- 
ations for the second sub-system with Zabs = 1.149092. A 
larger dispersion of At; in this case (peak-to-peak varia- 
tions ~ 300 m s~^ , except the 6th exposure) is a conse- 
quence of the lower accuracy of the line centering because 
of considerably smaller equivalent widths of the Fell lines. 
Table 7 lists the model parameters for the second Fe ii en- 
semble where rig = 2 (p = 29). 

Fig. 5 (^;abs= 1.150965) and Fig. 6 (2abs= 1.149092) 
show the synthetic profiles (smooth curves) with corre- 



sponding QSO data (normalized intensities) plotted by 
histograms. The positions of the sub-components are 
marked by vertical dotted lines. The corresponding X^j^ 
values are 0.9 {v = 2642) and 0.8 (i/ = 399) for the former 
and the latter sub-systems, respectively. 



4. Aa/a measurements 

The relativistic corrections of the Fe ii transition frequen- 
cies to the changes in a (the so-called g-coefficients) have 
been calculated by Dzuba et al. (2002). In our analysis we 
use the dimcnsionlcss sensitivity coefficients, Q = q/uiQ, 
listed in Table 5 (here = 1/Ao is the laboratory 
wavenumber). Then, the value of Aa/a can be estimated 



S. A. Levshakov et al.: Most precise single redshift bound to /^a/a 



9 



Table 6. SIDAM analysis: relative positions of the Fell lines (with respect to the adopted redshift of 1.150965), 6v, 
the line broadening velocities, h, the column densities, log A'', and the corresponding rms errors (given in parentheses) 
of the 13 components constituting the Fell absorption complex at ^abs = 1-151 





Positions of the main subcomponent in 


Fell 




different exposures, v (km s ) 




line (A) 


#1 #2 


#4 #5 


#6 #7 


1608 


U.0OD(U.U0Dj U.0ii(U.U80j 


0.9U((U.UOOj U.OOZ^U.Oooj 


U.dod(U.UodJ U.oiU^U.Uo/ j 


2344 


U.Oao(U.U/Dj U.DOO(U.U(Dj 


U.(9U(U.U//j U.oU4(U.U//j 


A OA/l/n ATTA n C/IC/'A A'7'7^ 

U.oU4(U.U/^ /^j U.040(U.U/ () 


2374 


U.OOz(^U.Uo/j U.DOD(U.Uo4j 


0.< 19(^0. 0o6j 0.(zy^U.Oooj 


A ctd'y ( r\ Aoc;\ n ^jaa/^a noc\ 

U.DD / (U.Uoo ) u.ouyt^u.uoo j 


2586 


U.498(U.Uf5Uj U.DiO^U.UoUj 


U.0/9(U.UoUj U. /o9(U.UoUj 


A Q1 ofc\ AQA^ n ccyiyi/'A n'7n\ 
U.olo(^U.UoU ) U.d44(^U.U / 9 ) 






Model parameters 




Component 


Sv, 


b, 


log AT, 


No. 


km s~^ 


km s ^ 


cm ^ 


1 


0.000 


2.55(0.13) 


13.243(0.049) 


2 


-39.23(0.10) 


3.71(0.08) 


12.701(0.080) 


3 


-30.02(0.67) 


4.07(0.76) 


12.790(0.079) 


4 


-24.96(0.19) 


2.31(0.09) 


13.458(0.075) 


5 


-20.56(0.40) 


2.01(0.33) 


12.859(0.099) 


6 


-14.95(0.17) 


2.33(0.17) 


13.245(0.081) 


7 


-3.85(0.37) 


5.05(0.15) 


13.268(0.111) 


8 


13.74(0.09) 


7.33(0.12) 


13.205(0.068) 


9 


28.46(0.37) 


4.69(0.49) 


12.893(0.111) 


10 


35.23(0.14) 


3.03(0.31) 


12.987(0.102) 


11 


44.04(0.20) 


4.35(0.25) 


13.287(0.078) 


12 


50.68(0.14) 


3.96(0.10) 


13.313(0.080) 


13 


58.12(1.22) 


18.53(0.93) 


12.913(0.054) 



Table 7. Same as Table 6 but for the Fell absorption complex at 2:abs = 1.149092 (the two component model) 



Positions of the main subcomponent in 
Fell different exposures, v (km s"'^ ) 

line (A) #1 #2 #4 #5 #6 #7 



1608 


-0.821(0.557) -0.988(0.572) 


-0.854(0.413) -0.695(0.376) 


-1.457(0.319) -1.010(0.226) 


2344 


-0.989(0.112) -0.867(0.096) 


-0.739(0.104) -0.978(0.160) 


-0.762(0.121) -0.914(0.152) 


2374 


-0.902(0.109) -0.713(0.245) 


-1.014(0.212) -0.595(0.198) 


-1.325(0.255) -1.165(0.119) 


2586 


-0.928(0.088) -1.017(0.069) 


-0.683(0.183) -1.068(0.167) 


-0.828(0.106) -1.245(0.106) 






Model parameters 




Component 






log AT, 


No. 


km 


km s ^ 


cm~^ 


1 


0.000 


1.51(0.01) 


12.501(0.002) 


2 


3.39(0.02) 


5.60(0.01) 


12.536(0.002) 



from a pair of lines with different sensitivity coefficients. 

In linear approximation (|Aq;/q;| <C 1), eq.(5) from LCMD 
can be re-written in the form: 

a 2c(Qi-Q2)' 

where the line positions Vi and V2 are taken from the same 
exposure. Here index '1' is assigned to the line A1608, 
while index '2' marks one of the other Fell lines (A2344, 
A2374, or A2586). The calculated Aa/a values are given 
in Table 8 along with their rms errors estimated by the 
standard method of error propagation. 



We now have everything we need to compute statis- 
tics: the mean value (Aa/o;) and its error. Wc note that 
for each individual exposure the uncertainty of Aa/a is 
dominated by the error of the line centering of the 'blue' 
Fell A1608 line since its strength is smaller as compared 
with the 'red' lines and the Fe ii A1608 spectra have sys- 
tematically lower S/N. Our approach reveals also a sys- 
tematic effect for all the exposures except the first one for 
the subsystem Zabs = 1.150965: the individual values of 
Aa/a have the same sign. At ^abs = 1.149092, the sign of 
Aa/a does not vary in the 1st and 6th exposures, whereas 
different sings are observed in the 4th, 5th, and 7th expo- 
sures. However, being combined the whole sample provides 
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Table 8. SIDAM analysis: Aa/a (ppm) calculated with eq.||2Jl and their drms errors (given in parentheses) 



Exp. 




Fe 11 pairs 




No. 


1608/2344 


1608/2374 


1608/2586 






= 1.150965 




1 


0.07(3.92) 


-0.73(3.45) 


1.03(3.27) 


2 


-6.05(3.90) 


-4.12(3.38) 


-2.97(3.26) 


4 


3.99(3.93) 


5.32(3.44) 


6.34(3.26) 


5 


-7.57(3.90) 


-4.17(3.37) 


-4.36(3.24) 


6 


-5.74(3.92) 


-0.88(3.41) 


-4.93(3.26) 


7 


5.61(3.95) 


5.68(3.42) 


7.39(3.27) 




•2abs 


= 1.149092 




1 


5.69(19.34) 


2.28(16.05) 


2.95(15.68) 


2 


-4.12(19.74) 


-7.77(17.59) 


0.83(16.02) 


4 


-3.93(14.51) 


4.51(13.14) 


-4.76(12.57) 


5 


9.64(13.91) 


-2.84(12.02) 


10.37(11.44) 


6 


-23.57(11.61) 


-3.73(11.54) 


-17.47(9.34) 


7 


-3.27(9.27) 


4.40(7.23) 


6.55(6.95) 



19 negative and 17 positive Aa/a values which random- 
ize these systematic shifts. Moreover, from inspecting the 
data listed in Table 8, one finds that two values —23.569 
and —17.466 from the 6th exposure (^abs — 1.149092) are 
sharply distinct from the others and exceed Scrims from the 
weighted mean. These measurements were rejected from 
the further analysis. Thus the final ensemble consists of 
the equal number {n — 17) of negative and positive Aa/a 
values. 

For the weights, Wi, we use a~^^ values from Table 8. 
The unbiased estimators of the weighted mean x = 
(Aa/a), and the variance of x are the quantities (Linnik 
1961): 



(4) 



\w\ 



and 

Var{x) 



\waa\ 



(5) 



(n — l)[w] ' 

where d — Xi — x are 'apparent errors', and square brack- 
ets denote summation, e.g., [w] — X]"=i '^iy i'^^] = 

En 
i^l WiX,. 

We find that the weighted mean of the ensemble of 
n = 34 Aa/a values and the accuracy of its determination 
are equal to (Aa/a) = -0.07 ± 0.84 ppm {la C.L.). 

5. Discussion and Conclusions 

We re-analyzed the profiles of the Fell lines associated 
with the sub-DLA system observed at Zabs = 1-15 in the 
spectrum of HE 0515-4414. Our main purpose was to 
achieve the highest possible accuracy in the line centering 
to set the most stringent constraint at a single redshift to 
the hypothetical variability of the fine-structure constant. 

This work was inspired by our recent results (LCMD) 
showing the influence of reduction errors and instabili- 
ties in the instrument on the relative radial velocity shifts 



in the merged QSO spectra. It was demonstrated that 
co-added spectra, minded to enhance the signal-to-noise 
ratio, may affect the shapes of narrow absorption lines be- 
cause each calibrated exposure has its own velocity offset. 

In the previous estimations of Aa/a from the co- 
added spectrum of HE 0515-4414 (QRL), we found that 
in spite of a good fit of the 8-componcnt model of the Fe ii 
profiles to the high S/N data, the normalized residuals 
showed an unexpected pattern in their radial velocity dis- 
tribution (see Fig. 7, upper panel). It was suggested that 
this pattern may be caused by the presence of unresolved 
narrow lines (or non-Gaussian line profiles) which may in- 
duce some kind of correlations between the residuals. We 
now recognize that additional correlations may be due to 
zero-point errors as well. 

If individual scientific exposures are treated separately 
and the number of the model components is increased 
up to 13, the pattern is vanished and the residuals be- 
come uncorrelated (see Fig. 7, lower panel). This demon- 
strates a clear advantage of the differential measurements. 
Moreover, the precision of the new estimate of (Aa/a) is 
improved by a factor of 2. From a general point of view, 
such an improvement would require 4 times longer total 
exposure time under the same observational conditions 
(i.e., 36h vs 9h, in our case). 

In this respect, the present work indicates that all steps 
in the data reduction procedure and in the long-term sta- 
bility of the instrument must be of a particular concern 
while dealing with sub-pixel positional measurements with 
echelle spectrographs. To make full and accurate utiliza- 
tion of the information derived from the observations a 
detailed knowledge of the instrumental characteristics is 
required. 

In the present analysis, we have reached an accuracy 
in the line centering which is comparable to the accuracy 
of the wavelength scale calibration, tTrms ~ 1 mA. From 
eq.© it follows that this value corresponds to the sys- 



S. A. Levshakov et al.: Most precise single redshift bound to /S.a/a 



11 



250 



- 



-250 



2. 500 



T I I r 




Exposures: 
1 



z.,, = 1.150965 
J I I I 



250 



-250 



-500 



-750 



T 1 1 r 




Exposures: 
1 



■ I 




Zabs = 1 


149092 








. 1608 


2344 


2374 


2382 


2586 


2600 


A 


- 1.97 
1 


2.43 
1 


1.87 
1 


2.88 
1 


2.25 
1 


2.79 
1 


log fA 
1 


19b 


Br 


9r 


lOr 


19r 


20r 





echelle order 

Fig. 4. Same as Fig. 3 but for the Fell lines from the 
Zahs = 1.150965 and 2;abs= 1.149092 systems. The corre- 
sponding models are superpositions of Ug = 13 and rig = 2 
Voigt profiles. The upper panel illustrates the effect of 
hidden (unresolved) components: the peak-to-peak varia- 
tions of Aw are decreasing with increasing number of sub- 
components. Smaller equivalent widths of the Fell lines 
at Zahs= 1.149092 lead to a larger dispersion of the At; 
values shown in the lower panel 



tcmatic error of the individual determination of Aa/a 
of CTgys ~ 4 ppm (calculated from a pair of Fell lines 
with |AQ| = 0.06). This adds a systematic error of about 

1 ppm to the total error budget of the weighted mean, 

^lot = ^(A«/a> + <s - 1-19 PPm^ 

It should be emphasized that the level of {Aa/a) ~ 
(—5.7 ± 1.1) ppm found by Murphy ct al. (2004) was es- 
timated from a large sample of 143 absorption systems 
ranging from z = 0.2 to 4.2 which were observed with the 
HIRES/Keck spectrograph (here the error of the mean 
includes only statistical uncertainties). Since for the in- 
dividual system at z = 1.15 the accuracy of (Aa/a) is 
better than 1 ppm, we are now in the position of arguing 
if the Keck ensemble average is biased. 

We notice that the level of Aa/a expected at 2; = 1.15 
from the damped-oscillatory model by Fujii (2005) is 

2 ppm (see Fig. 8). We do not see this value in our analysis. 
However, the total error of (Aa/a) is not small enough to 
verify or reject Fujii's model. With the probability level 
of 0.05 we cannot take the observed difference between 
—0.07 ppm and 2 ppm as significant according to the t- 
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Fig. 7. The normalized residuals, (JT*^"' — T"'"') /a, for the 
Fell profiles (AA1608, 2344, 2374, and 2586 A) calculated 
from the 8-component model in QRL {upper panel) and 
from the 13-component model in the present work {lower 
panel). The zero radial velocity is fixed at 2; = 1.150965 

test. To probe the oscillatory behavior of a, very accurate 
measurements of Aa/a at higher redshifts (where the am- 
plitude of Aa/a is expected to be ±5 ppm) are required. 

As a conclusion, it is worthwhile to note that the 
achieved accuracy of Aa/a is unique for the standard 
UVES configuration and that further improvement at the 
sub-ppm level can be attained only with increasing spec- 
tral resolution and stabilizing instrumental performance 
such as, for instance, a fiber link producing a stable illu- 
mination at the entrance of the spectrograph and allowing 
continuous simultaneous comparison spectrum. 
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